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Abstract 

Within a dynamical coupled-channels model that has already been fixed by analyzing the data 
of the irN —* ttN and *yN — > ttN reactions, we present the predicted double pion photoproduction 
cross sections up to the second resonance region, W < 1.7 GeV. The roles played by the different 
mechanisms within our model in determining both the single and double pion photoproduction 
reactions are analyzed, focusing on the effects attributable to the direct "yN — ► ttttN mechanism, 
the interplay between the resonant and nonresonant amplitudes, and the coupled-channels effects. 
The model parameters that can be determined most effectively in the combined studies of both 
the single and double pion photoproduction data are identified for future studies. 

PACS numbers: 13.75. Gx, 13.60.Le, 14.20.Gk 
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I. INTRODUCTION 



The spectrum and structure of low-lying nucleon and A resonances (collectively referred 
as iV*) are primordial information for any understanding of the nonperturbative QCD do- 
main. Consequently, a great effort has been made at the Excited Baryon Analysis Center 
(EBAC) during the past few years to extract the properties of N* from the world data on 
ttN -> tcN and 7 A -> ttN data Q. 

It is well acknowledged nowadays that a proper extraction and further interpretation of 
iV* properties require the construction of reaction models that maintain the unitarity of most 
relevant channels and can correlate the vast amount of data for both the single and double 
meson production reactions. Among the existing theoretical approaches, the one taken 
at EBAC tries to encompass the aforementioned by considering the interactions among the 
7JV, ttN, r]N, and itttN channels within a multichannel,multiresonance framework J2[. After 
constraining the hadronic part of the model by fitting [3J the ttN — > irN scattering data, we 
have performed our first studies of single pion photoproduction [4j and electroproduction 
reactions 

As discussed in our previous works, the hadronic part of the model was constrained 
mostly using tcN —>■ ttN experimental data. This means that the couplings of the N* to 
the 7rA, pN and aN channels, which are the quasi- two-body channels of the tittN, are 
necessarily not well constrained in the current version of the model. To this extent, double 
pion photoproduction reactions are important for understanding the way iV* couple to the 
7nrN channel, and thus to refine our global dynamical coupled-channels framework. In 
Ref. 0], we carried out such a study for nN — > rniN reactions with the predicted cross 
sections in reasonable agreement with the available data. In this work, we extend that work 
to investigate double pionphotoproduction reactions by comparing our predictions with the 
total cross sections data 0, @, 0, 10, [ll|, 12, 13] and invariant mass distributions 



We first present the predictions of our model for the double pion photoproduction reactions 
up to W = 1.7 GeV. We then analyze how the discrepancies with the data are sensitive 
to which of the electromagnetic parameters of the model, as a step toward performing the 
combined fits of the world data on irN, 7 A — > ttN, tttcN reactions. 

Most of the previous investigations of the double pion photoproduction reactions em- 
ployed the tree-diagram models [H, 0, 17, T^|, emphasized the roles of certain resonances 



on specific double pion photoproduction reactions, or focused on the very near threshold 
region using chiral perturbation theory (l9l . 2oT ]. In our approach, we do not make such 



simplifications. We perform the full coupled-channels calculations and include all channels 
and N* states determined in Refs. [!, HJ. 

The basic formulas used in this work are presented in Sec. [Til In Sec. IIHI we present the 
predictions of the current model and analyze the contributions from the direct 'jN — > inrN 
mechanism and the transitions from 7 A to the unstable 7rA, aN and pN states. In Sec. IIVI 
we scrutinize the contribution of each of the 7iV — > N* helicity amplitudes on both single 
pion and double pion photoproduction reactions. A summary and some conclusions are 
given in Sec. [V] 



II. BASIC FORMULAS 

Within the EBAC dynamical coupled-channels (EBAC-DCC) model, the 7A — > tyttN 
amplitude consists of four pieces (see Fig. HD: 
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FIG. 1: Graphical representations of T 7 jv )7 r7rjV of Eqs. (HJ-©. 

T lN ,TmN(E) = T^ nnN (E) + T^ nnN (E) + T? N}7rnN (E) + T°$ )WirN (E), (1) 

with 

T^n,ttttn(E) = v jn,ttttN + ^2 T 7 n,mb(E)Gmb(E)vmb,7vwn , (2) 



MB 



with 



T*n,itwn(E) = T 1 n,tta(E)G 7T a(E)Ta^ttN, (3) 
T in,ttkn{E) = T 7N)PN (E)G pN (E)h p ^- K7r , (4) 

r EyN,TnrN(E) = T lNtCrN (E)G a N(E)h (7 ^ 7r7T . (5) 

Here Ta^ttAt, h p ^ n7T , and h a ^ wn describe the A — > nN, p — ► 7T7r, and a — > 7r7r decays, re- 
spectively; Gmb{E) (MB = irN,rjN,7rA, pN,aN) are the meson-baryon Green's functions. 
v-yN^-wN represents the direct — > nnN transition potentials illustrated in Fig. [2j The 
processes described by v^n^-kN are not contained in the T^ B n7rN , and thus there is no double 
counting. 

The 'jN — > MB transition amplitudes can be divided into the so-called nonresonant and 
resonant amplitudes (suppressing angular momentum, isospin, and momentum indices), 

T y N,MB(E) = t 7 N,MB{E) + t^ NMB (E), (6) 

\ LA Nr.'' l \.-' l LJ^iZ 
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FIG. 2: Diagrams considered for v^n^-kN- 
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with 

t^N,MB{E) = V lN ^[ B + V 1 N t M'B'GMiB'{E)t M i B ^ MB {E). (7) 

M'B' 

and 

tyN,MB( E ) = Yl r jN ^ N *(E)[D(E)] itj T N *->MB(E) . (8) 

TV* ,7V* 

In Eq. ((Tj), v.yN,MB represents the 7iV — > M5 transition potential derived from tree diagrams 
of a set of phenomenological Lagrangians describing the interactions among 7, tt, rj, p, cu, 
cr, N, and A(1232) fields. The details are given explicitly in Appendix F of Ref. [2J. The 
dressed 'yN — > N* vertex function appearing in Eq. ([HD is defined by 

FjN-*n* (E) = r 7 jv-,iv* + X! v ^n,m'B'Gm'B'(E)Tm'B'^n*(E), (9) 

M'B' 

where T^n^n* denotes the bare jN — > N* vertex within the EBAC-DCC model and is 
parametrized as 

r^Ajvte) = 7^^\/^^^[v / 2^ A a]^,(a 7 -a jv ), (10) 

where q R is defined by the N* mass M N * = qn + E N (q R ). 

Within our model, the meson-baryon Green function Gmb, the hadronic nonresonant 
amplitude t MB ^ M / B >, the dressed N* propagator D(E), and the dressed iV* — > MB vertex 
function T N *^ MB are purely hadronic processes. We take these hadronic pieces from the 
model constructed from analyzing the data of ttN — > ttN scattering [^], and keep them fixed 
throughout this paper. 

The calculation of the terms T^^ nN with MB = nA, pN, aN, defined by Eqs. ([3])- 
([5]), is straightforward. However, the calculation of the second term of T^§^ nnN , defined by 
Eq. ([2]), is much more complex. To simplify the calculation, we employ the same prescription 
as in the calculation of the ttN — > ttttN reactions 0]. This is based on the observation that 
the processes illustrated in Figs. [2](a)-12](d) can be written as 

V tN\1-kN ~ v -yN,TTNG 7T N(E)h N ^ 7rN , (11) 

where v^fl^. N is the sum of the all processes illustrated in Figs. |2^a)-|2](d), v^n^n is the 
two-body 7iV — > ttN potential, and h^^N is the iV — > ttN vertex function. Taking 
account of only a part of vmb^-kN that can be approximately expressed as vmb,itkN ~ 
VMB,TrNG w N(E)h N ^ 7rN , Eq. ([2]) can be written as 

T^N,irnN(E) ~ V ^N,TrnN + [ v ^N,-nN + T^n^mb(E) G m b{E)vm B ,n n]G n N {E) h n N 

MB 

= V ^N-kttN + T^N,TrNG 7T N(E)h N ^ 7rN . 

(12) 

Here in the last step we have used the relation T iN ^n = v iN ^n + Y,mb T^n,mbGmbVmb,ttN- 
We use Eq. (1121) which can be calculated with all parameters taken from our previous analysis 
of irN, 7iV — > ttN reactions. 

The formulas for calculating total cross sections and invariant mass distributions from 
our amplitudes can be found in Ref. [6| and are not shown here. 
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FIG. 3: Near threshold behavior of the total cross section for ^yp — ► ttttN: (a) 71? — ► ir + ir~p, (b) 
7P — > 7r°7r°p, and (c) 7p — > 7r + 7r°n. The red solid curve is the full results predicted from our current 
model, and the blue dashed curves are the results without the n7vN contribution. The data are 
taken from Refs. 0, fi 0, M, EJ Q Q • 



III. ANALYSIS OF THE DIRECT REACTION MECHANISMS AND THE 
COUPLED-CHANNELS EFFECT 

With the parameters determined from our previous analysis of itN, — > ttN reac- 
tions 0, 0], the results presented in this section are pure predictions within the current 
model developed in EBAC. We first present our results of the double pion photoproduction 
reactions, and then examine how the reactions mechanisms within our model determine the 
cross sections. 

In Fig. [3X we find that our current model (red solid curve) has a good agreement with 
the 7iV — > ttttN total cross sections in the energy region up to W — 1.4 GeV. We observe 
that the direct T^ nwN amplitude can greatly improve the model to reproduce the near 
threshold behavior of the jN — > ttttN total cross section data. Its effects in higher W are 
shown in Figs. SJa), SJb), and|H(c). The red solid curves are the predictions from our full 
calculations and the blue dashed curves are from turning off the term T^^ n7rN (the bands in 
the figure will be explained later in this paper). We see that the effect of T^ n7rN is sizable 
on 7£> — > 7T + 7T~p [Fig. H](a)] and ^p — > 7r°7r°p [Fig. H](b)], and negligible on 7^ — > 7i + 7T°n 
[Fig. BJc)]. It is clear that its inclusion does not change the energy dependence of the total 
cross sections for any of the considered ^yN — > ttttN reactions. 

Although the threshold behavior is in general well reproduced as can be seen in Fig. [3l 
our predictions at higher W shown in Fig. H] clearly overestimate the experimental data 
above W = 1.4 GeV in both 'jp — > tt + tt~p and 7p — > 7i°n°p reactions, while the results of 
7£> — > 7i + 7i°n are good up to W — 1.5 GeV. However, our current model reproduces the 
7iV — > ttN reactions quite well in the considered energy region, as seen in the right panels 
of Fig. HI This fact indicates that there exist reaction processes which have significant effect 
on the observables of 7iV — > tcttN, but not of 7iV — > tcN. 

To get some insights into our disagreement with the data and to guide our future combined 
analysis of all nN, 'jN — > ttN, ttttN reactions, we examine which mechanisms are most 
relevant to our calculations in this energy region. We first examine the contributions of 
each process appearing in Eqs. (J2])-(poD- The results from the full amplitude are shown in the 
top row of Fig. EJ T*$ >7nrN (black solid), T°$ >7nrN (red dashed), T$^ nN (green dotted), and 
T^NnirN (blue dash-dotted). The figures in the left, middle, and right columns are of the 
7P — > 7t + tt~p, '-/p — > 7T°7T°p, and 7p — > 7T + 7T°n total cross sections, respectively. 

We also show in the middle (bottom) row of Fig. [5] the results for which the full two- 
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FIG. 4: Total cross sections of the double and single pion photoproduction reactions up to W = 1.7 
GeV: (a) ^p — > ir + TT~p, (b) ^p — > tt tt p, (c) *yp — ► 7r + 7r°n, (d) 7p — ► 7r°p, and (e) ^p — ► TT + n. The 
red solid curve is the full result predicted from our current model, and the blue dashed curve in (a)- 
(c) is the result without T^§^ nN contribution. The band is generated by allowing a 25% variation 
in the value of the 7riVA coupling constant g^NA used in the electromagnetic amplitudes. The data 
of the double and single pion photoproduction reactions are taken from Refs. [?| la.|glEo. 
and Refs. [2lfj , respectively. 
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body amplitude T lN>MB in Eqs. ©-(JSD an d (TT21) is replaced with its resonant (nonresonant) 
part T 7 tv,mb - * t^N,MB {T~/N,mb ~~ > ^n,mb)- Thus we can examine the relative importance 
between different mechanisms in resonant t^ N MB and nonresonant t lNiMB amplitudes sep- 
arately. Note that the curves describing the 7iV — > aN ( / ~fN — > pN) process are not seen 
in the jp — ► 7T + 7i°n (jp —>■ 7r 7r°p) total cross sections because the corresponding terms 
do not contribute because of isospin selection rules. In Figs. Eta)-|S(c), we clearly see that 
the full 7iV — > 7rA — > ttttN processes (black solid curves) have the largest contribution 
compared to the other processes. By comparing Figs. 0(a) and 0(g), we further find that 
the large discrepancy with the 7^ — > tt + 7t~p data is attributable mainly to the nonresonant 
7iV — > 7i A — > ttttN amplitude. The dominance of the nonresonant 7iV — > 7rA — > 7T7tA^ in all 
three 7iV — > ttttN reactions can also be seen in the bottom panels of Fig. [5j 

Most of the nonresonant 7iV — > 7rA transition matrix elements considered in our model 
depend on the 7riVA coupling constant g n NA (see Ref. for the details). We thus examine 
how our predictions are sensitive to this coupling strength. This is illustrated in Fig. [3] where 
we have presented bands, which are generated by varying g W NA included in the 7iV — > nA 
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FIG. 5: Contributions of each reaction process described in Eqs. d2J)- dSJ) to the total cross sections. 
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transition matrix elements by ±25%. Clearly such changes in g^NA have a great influence 
on yp — > 7r + 7r~p (top) and yp — > n + n°n (bottom), and less of an influence on yp — > Tc°7i°p 
(middle). Within our dynamical coupled-channels model, the 7iV — > 7rA process also enters 
in the single pion photoproduction reactions as a consequence of the unitarity, and thus its 
change consistently affects the single pion photoproduction observables, too. As can be seen 
in the right panels of Fig HI its importance turns out to be very minor in the 7iV — > ttN 
total cross sections. The bands from varying g n NA in 7^ - * ttA by ±25% are not visible. 
From this observation, in the remainder of this paper we will use a 20% smaller value for 
the g n NA appearing in the electromagnetic potentials. The value turns out to be very close 
to that of the quark model. 

In Figs. [MHl we show the predicted invariant mass distributions of yp — ► 7T + 7T~p, 
yp —>■ TT°ii p, and yp — > 7r + 7r°n, respectively. To compare with the shapes of the data, 
the overall magnitudes of our predictions (red solid curves) are normalized to have the same 
integrated values of the data. We can see that the shapes of the predicted nN invariant 
mass distributions are in reasonable agreement with the data for all cases considered, while 
deviations are seen in several tttt invariant mass distributions (right panels of Figs. [ME])- 

This is found to be attributable to the fact that the irN distributions are dominated by 
the A(1232) in the 7iV — > 7rA(1232) — > ttttN process, while the irir distributions involve the 
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FIG. 6: Invariant mass distributions of ^p — > ir + Tr~p at W = 1420,1520 MeV: (left) (tt~p); 
(middle) (vr + p); (right) (ir + ir~). The red solid curve is the full result, and the blue dashed curve is 
the phase space distribution. The magnitude of both curves is normalized to the data. The data 
are taken from Ref. 1141. 



interferences among all of the *yN — > irA, pN, aN — » ttttN amplitudes. The results of the 
7T7T invariant mass distributions have provided useful information for improving our current 
model. In particular, the deviations from the data in the 7t + tt° distributions of the TT + -K°n 
channel at high invariant mass (right panels of FigJH]) suggest that the parameters associated 
with the pN channels will need to be modified. 

The most relevant novelty of the present study is the use of a dynamical coupled-channels 
model. In Fig. M, we show the coupled-channels effects associated with the electromagnetic 
interactions on the jN — > tittN total cross sections, which is demonstrated here for the first 
time in the investigations of double pion photoproduction reactions. The red solid curves 
are our full results. The green dotted curves are the results in which only the diagonal part 
(M'B' = MB) is taken in the M'B' summation of Eqs. ([7]) and ([9]), and the blue dashed 
curves are obtained by further setting t 7 v,A/s — > v^n,mb and T^n^n* — > ^-yN^N*- These 
correspond to examining the coupled-channels effect associated with the electromagnetic 
interactions. (Note again that the pure hadronic part of the amplitudes is fixed with the 
model determined in Ref. 0] throughout this paper.) In the considered energy region up to 
W = 1.7 GeV, we find that the blue dashed and green dotted curves almost overlap with 
each other but both of them are quite different from our full results (red solid curves). This 
suggests that the structure in the 7p — > n+Ti^p, n 7T p total cross sections is attributable 
mostly to the couplings between reaction channels. 

Before closing this section, we comment on the recent measurements of the polarization 



observables. It was shown in Refs. [22|, |23| that existing reaction models have significant 



discrepancies in the beam-helicity asymmetry measured at CLAS 22j and more recently at 
MAMI [23| • We have observed that our current model also produces similar discrepancies to 



that of other works shown in Refs. [22l. |23| . These results indicate that the polarization ob- 
servables will provide critical information on constraining reaction models and understanding 
the iV* states. 
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FIG. 7: Invariant mass distributions of -yp -> 7r°7r°p at = 1420, 1520 MeV: (left) (vr°p); (right) 
(7r°7r°). The red solid curve is the full result, and the blue dashed curve is the phase space 
distribution. The magnitude of both curves is normalized to the data. The data are taken from 
Ref. 0]. The energy bins of the data are 20-30 MeV around the central W shown in the panels. 
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FIG. 8: Invariant mass distributions of 7p 7r + 7r°n at W = 1420, 1520 MeV: (left) (vr°n); (middle) 
(-7r + n); (right) (7r + 7r°). The red solid curve is the full result, and the blue dashed curve is the phase 
space distributions. The magnitude of both curves is normalized to the data. The data are taken 
from Ref. [lOj]. The energy bins of the data are 20-30 MeV around the central W shown in the 
panels. 



IV. EFFECTS OF RESONANCES 

The bare helicity amplitudes, defined in Eq. (IT0I) . are free parameters in our framework. 
They quantify the photoexcitation of the core N* states and, together with their dressed 
counterparts, are to be interpreted by means of microscopic models (e.g., quark models, 
lattice QCD calculations). Although A J X are taken to be real numbers, the dressed helicity 
amplitudes, which have in general a sizable contribution from the second term in Eq. (Q, are 
complex numbers. This second term contains the meson-cloud contribution to the 'jNN* 
vertex, which is to a large extent fixed from the strong interaction sector. 

In this section we present the effect on the single and double pion photoproduction 
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FIG. 9: Coupled-channels effects associated with electromagnetic interactions. The red solid 
curve is full results, the green dotted curve corresponds to taking only the diagonal element in the 
M'B' summation in Eqs. ([7]) and ([9]), and the blue dashed curve is obtained by further makinga 
replacement of iyjv.MB — ► v^n^mb and F^n-^N* IVv^tv*- The data are taken from Refs. 0,M 
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observables of variations on the bare helicity amplitudes, which affect directly the dressed 
ones, see Eq. (Q. This will be done by presenting results computed by varying the initial 
value of the bare helicity amplitudes listed in Tabled! by ±50%. The results are presented 
as bands in the figures, these bands are generated by filling the region enclosed by curves 
from two calculations using 0.5 x A,/ 2 and 1.5 x Aj/ 2 . 

Before proceeding to showing our results, we comment on the bare helicity amplitudes 
presented in Table [H Those values are not exactly the same as those of our previous 7iV — > 
ttN analysis [3] • There we did not provide any measure of the uncertainty in the bare helicity 
amplitudes that resulted from fitting the photoproduction data. In the current paper we have 
varied the binning of the data and thus some of the less constrained helicities resulting from 
the fit are varied. In the following we will quantify the effect of such variations, providing a 
clear indicator of the dependence of our results for both single and double pion production 
on the helicity amplitudes. 



TABLE I: The bare 7./V — > N* helicity amplitudes determined from x 2 -fits to the jN — > 7riV. The 
asterisks in the second (third) column mark the iV* states in which jN transition process is found 



to be relevant to the single (double) pion photoproduction reactions up to W = 1.7 GeV. 
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FIG. 10: The panels (a)-(e) depict the total cross section, differential cross sections and photon 
asymmetry for ~/p — > Tr°p, and (f)-(j) show total cross section, differential cross sections and photon 
asymmetry for ^p — > 7r + n. Each band is obtained by allowing a 50% variation of the helicity 
amplitudes for the Ai/ 2 of 5n(1535) (solid blue) and A\i 2 of 53i(1620) (oblique-lined red) listed 
in Table HI The data are taken from Ref . [2lf] . 



A. S-wave iV*s 



We start the comparison with the Sii(1535) and S , 3i(1620). In Fig. [10]we show the effect 
of varying their helicity amplitudes on the single pion photoproduction data. The sample 
data we consider are the total cross sections (left panels) for 7p — > 7r°p and 7p — > ix + n and 
differential cross sections and polarization data in the A(1232) region (middle panels) and 
in the W = 1500 MeV region (right panels). 

First we note that the ±50% change in helicity amplitudes for the S , n(1535) resonance 
plays an important role in building the peak near the 1500 MeV region for both 7^ — ■> 
7r + n, 7r°p total cross sections[see Figs.fTOTa) andfTDTf)] and correspondingly in the differential 
cross section near the 1500 MeV region [see Figs.fTOTc) and fTDT h)]. The S31 gives a prominent 
contribution in the whole energy region above the A (1232) region, as indicated by the 
oblique-lined bands. The S31 also affects the forward peaking of the 7^ — > 7r + n differential 
cross section data around W = 1500 MeV [see Fig. fTDT h)]. Their influence on the photon 
asymmetry S is sizable and qualitatively similar for both resonances, being negligible in the 
A(1232) region. 

Now we turn to the double pion photoproduction reactions, see left column of Fig. [Til 
First, as expected, and the same occurs for all resonances considered, the helicities have no 
influence on the near threshold behavior. Second, both S-wave resonances play a relevant 
role for the considered reactions. Modifying the Ai/ 2 of the 5*11(1535), the total cross sections 
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FIG. 11: Total cross sections, [panels (a), (d), and (g)] 7£> — > -K + ir~p, [panels (b), (e), and (h)] 
7P — > 7r 7r p, and [panels (c), (f), and (i)] jp — > ir + ir°n. The different bands are generated by 
allowing ±50% variations of the helicity amplitudes listed in Table [H The data are taken from 

Refs. B H El EJ EE Q • 

for 7p — > tt°tt°p and 7p — > 7r + 7r°n can vary up to 20%, although there is no qualitative change 
in the energy dependence of the total cross sections [see Figs. ITTTb) and fTPT c)]. The S31 case 
is similar, but actually affects all the reactions. A smaller value of the S31 helicity amplitude 
is suggested by these results. Within our model, none of the peaks seen in the total cross 
section data can be ascribed solely to S-wave resonances. 

B. P-wave N*s 

The helicity amplitudes of the A(1232) resonance are essentially fixed by analyzing data 
near its nominal mass, as has long been known. In Fig. [J2] we fully confirm this. The effect 
of variations on both Au% and A 3 / 2 of the A(1232) is well localized around its peak but 
reaches up to 300 MeV above it in the 7p — > ir°p reaction [see Figs. IT2T a)-[T2Te)]. This can 
also be seen in their influence on the photon asymmetry at W — 1480 MeV. The A3/2 mostly 
affects the perpendicular angles, while the A1/2 affects the forward and backward angles. 
The Roper resonance plays a minor role, with no sizable trace in the observables. 

In the double pion photoproduction case, however (see middle column of Fig. ITT]) , the jN 
transition processes of both P33 and Pu play almost no role in the entire considered region. 
Let us point out that we refer here to the influence of the P33 as an s-channel exchange, the 
importance of the A in this reaction is of course large, as pointed out in Sec. IHfl where we 
show that most of the reaction flows through the irA channel. 

C. D and F-wave iV*s 

Let us first study the influence of the helicity amplitudes on the single pion photoproduc- 
tion. The -D13 is responsible for part of the second peak near 1500 MeV in the 7p — > 7i + n 
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FIG. 12: The panels (a)-(e) depict the total cross section, differential cross sections and photon 
asymmetry for ~/p — > Tr°p, and (f)-(j) show total cross section, differential cross sections and photon 
asymmetry for ^p — > 7r + n. Each band is obtained by allowing a 50% variation of the helicity 
amplitudes for the A 3/ / 2 of ^33(1232) (solid brown), Ai/ 2 of ^33(1232) (oblique-lined red), and Ai/ 2 
of Pn(1440) (solid black) listed in Tabled The data are taken from Ref. [211 ] . 



total cross sections [see Fig. [TBTf)]. The F 15 (1680) contributes to the third peak in both 
total cross sections. In the middle and right panels of Fig. [13l we see that none of the ±50% 
changes of -D13, and -D33 helicity amplitudes affect much the £ and da/dfl observables. 

D wave resonances have long been advocated as being responsible for most of the structure 
observed in the total cross sections for ^p ttttN. The first peak in the total cross sections 
has been explained in tree level calculations thanks to the -Di3(1520) (lol . [lrl 18| and to 
interferences with the D 33 (1700) [13]. In our coupled-channels model we confirm the very 
important role played by the _Di3(1520), which builds up a large fraction of the first peak in 
the 7p — > ir n°p reaction [see the right panels in Fig. [11] . On the other hand its effect is also 
sizable on the 7p — > -k + -k^p total cross section, producing an overprediction of this observable 



in our model. As in the tree-diagram models of Refs. [151 . 1 161 . |18| the peak structure in this 
reaction is always much more pronounced in the models than in the experimental data. 
Effects of the _D 33 are sizable only on the jp — > 7r + 7r°n, similar to what was reported in 
Ref. ' 



17], but they do not produce a peak structure as the experimental data show. 



V. SUMMARY AND CONCLUSIONS 



Within the dynamical coupled-channels model constructed from analyzing the single pion 
production reactions [1, 0] , we have investigated the total cross sections and the invariant 
mass distributions for the double pion photoproduction reactions off the proton in the energy 
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FIG. 13: The panels (a)-(e) depict the total cross section, differential cross sections and photon 
asymmetry for ~/p — ► Tr°p, and (f)-(j) show total cross section, differential cross sections and photon 
asymmetry for ^yp — > 7r + n. Each band is obtained by allowing a 50% variation of the helicity 
amplitudes for the A s / 2 of Z?i3(1520) (solid brown), A s / 2 of Fi5(1680) (oblique-lined red), and 



A 3 / 2 of 1)33(1700) (solid black) listed in Table[IJ The data are taken from Ref. 21]. 



region up to W = 1.7 GeV. In the low-energy region up to W = 1.4 GeV, our results agree 
well with the total cross sections data, in which the direct process T^§ n7rN plays a crucial 
role for the reproduction of the data. Above W = 1.4 GeV, our current model starts 
to overestimate the data for yp — ► tt + 7t~p and yp — > 7r 7r°p. We have found that the 
jN — > ttA process is most relevant for the 7^ — > rnrN reactions and is a major origin of the 
overestimation in the 7p — > 7t + tt~p total cross section. Our model reproduces well the shapes 
of the invariant mass distributions data except for several 7T7T invariant mass distributions 
of 7p — > 7T°Tr°p and jp — > 7i + 7T°n. We expect that this deviation provides useful information 
to improve our current model. Also, we have demonstrated the coupled-channels effects on 
the double pion photoproduction case, which is of similar size to the 7riV — ► tttiN case. 

It is noted that our current model describes the single pion photoproduction observables 
in the same energy region quite well. We thus have examined the origins of our disagreements 
with the data by considering both the single and the double photoproduction reactions. We 
have found that the 7riVA coupling constant g n NA in the — > MB transition matrix 
element plays an important role. If we reduce its strength determined in Ref. [3| by 25 % to 
a value close to the quark model value, the magnitude of the ■yp — > n + n^p total cross section 
is drastically reduced, while the corresponding changes in the single pion photoproduction 
observables are negligible. This finding indicates that a smaller value of g W NA w iU be needed 
in a combined analysis of the world data of 7riV, yN — > tcN, tcttN reactions. 

We have also investigated the sensitivity of each 7^ — > N* process to the yN — > ttN 
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and 7 A -> rniN reactions. The 7 A -> S u (1535), 7 A -> S 3 i(1620) and 7 A -> L> 13 (1520) 
processes are found to have significant influence on both the single and the double pion 
photoproduction observables. In particular, 7 A — > £)i3(1520) will be key to fixing the 
overestimation at the first peak of 7 A — > -K°7i°p around W = 1.5 GeV. As for the P wave 
resonances, the 7A — > A(1232) process is critical for describing the 7A — > ttN observables 
up to W = 1.5 GeV, while it plays almost no role for the total cross sections and invariant 
mass distributions of 7 A — > ttttN reactions. The 7 A — > A* (1440) process just has a 
negligible contribution to the 7 A — > ttN, ttttN observables considered in this paper. This 
result for the A* (1440) is consistent with the recent analysis in Ref. (24J. The A* states that 
are found to be important in determining the single and double photoproduction reactions 
are indicated in the second and third columns of Table [B 

The results in this paper show clearly that in general the analysis of the single pion 
production reactions is not enough to pin down the amplitudes associated with the electro- 
magnetic interactions. To extract the reliable information on the A* states below W = 2 
GeV, at least one needs to perform simultaneous analysis of the single and double pion 
production reactions. Currently, this is one of the main efforts at EBAC. 
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